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We explore the optical properties of periodic layered media containing left-handed metamaterials. 
This study is based on several analogies between the propagation of light in metamaterials and charge 
transport in graphene. We derive the conditions when these two problems become equivalent, i.e., 
the equations and the boundary conditions when the corresponding wave functions coincide. We 
show that the photonic band-gap structure of a periodic system built of alternating left- and right- 
handed dielectrics contains conical singularities similar to the Dirac points in the energy spectrum 
of charged quasiparticles in graphene. Such singularities in the zone structure of the infinite systems 
give rise to rather unusual properties of light transport in finite samples. In an insightful numerical 
experiment (the propagation of a Gaussian beam through a mixed stack of normal and meta- 
dielectrics) we simultaneously demonstrate four Dirac point-induced anomalies: (i) diffusion-like 
decay of the intensity at forbidden frequencies; (ii) focusing and defocussing of the beam; (iii) 
absence of the transverse shift of the beam; and (iv) a spatial analog of the Zitterbewegung effect. 
All of these phenomena take place in media with non-zero average refractive index, and can be tuned 
by changing either the geometrical and electromagnetic parameters of the sample, the frequency and 
the polarization of light, or the frequency and the polarization of light. 



I. INTRODUCTION 



Similarities between Maxwell and Dirac 
equations 



Highly unusual properties of monolayers of graphite 
(graphene) and of optical media with negative refractive 
indices (left-handed metamaterials) had been indepen- 
dently predicted and studied theoretically a long time 
agnli^. At that time, however, these predictions were 
perceived as rather intriguing but unrealistic exotica, and 
remained unnoticed for about a half-century, until quite 
recently (and nearly simultaneously) they were embod- 
ied in real materials. This immediately triggered an ex- 
plosion of interest and activities, in metamaterials and 
graphene, both in solid states physics and optics. Re- 
searchers also realized that the most unusual proper- 
ties of electron transport in graphene were also pecu- 
liar to the propagation of light in dielectric systems with 
metamaterials. Mathematically, this is because, under 
some (rather general) conditions, the Maxwell equations 
for electromagnetic waves in an inhomogeneous dielectric 
medium can be reduced to the Dirac equations for charge 
carriers in graphene subjected to an external electric po- 
tential. 



The history of recasting Maxwell equations in alterna- 
tive, more compact, spinor forms goes back to the be- 
ginning of the past century, and is still in progress (for a 
comprehensive historical overview see Rcfs.landl, with 
recent examples in Refs.[H-[l). Therefore it is not surpris- 
ing that the similarity between Maxwell and Dirac equa- 
tions has long been noticed (according tc£, Majorana dis- 
cussed it already in 1930). In the general case of inhomo- 
geneous media, it was explicitly demonstrated in Ref. [H 
the close analogy between: (i) the quantum-mechanical 
form of the equations for the Reimann-Silberstain vector 
fields (linear combinations of the electromagnetic vectors 
D and B), and (ii) the Dirac equation, written in the 
chiral representation of the Dirac matrices was demon- 
strated explicitly in the general case of inhomogeneous 
media. 

Recently, as graphene became increasingly more popu- 
lar in solid state physics, the mathematically-established 
similarity of Maxwell and Dirac equations took on a new 
physical significance. Inspired by the very unusual pre- 
dictions and discoveries made in graphene, researches in 
optics started endeavors to reproduce the unique trans- 
port properties of graphene in specifically-designed di- 
electric structures. An additional incentive to these ef- 
forts came from the fact that, while the elementary build- 



ing blocks of graphene are fixed, modern micro- and 
nanotechnologies enable manufacturing periodic dielec- 
tric samples with a variety of types and sizes of unit 
cells. Moreover, the electrodynamic parameters of pho- 
tonic crystals can, in principle, be controlled by external 
fields, providing unique opportunities to study condensed 
matter phenomena in optical ways; for example, by open- 
ing a gap between Dirac cones, as well as breaking and 
restoring space-inversion and time-reversal symmetries 5 . 
Furthermore, rather simple electrodynamical analogies 
furnish physical insights into properties and applications 
of graphene such as: the Klein phenomenon^, breaking 
the valley degeneracy-^, graphene quantum dot a 11 ' 12 , the 
electronic Goos-Hanchen shift^, derealization in one- 
dimensional disordered systems^, etc. The implemen- 
tation of the analogy between Dirac electrons and light 
could be rewarding for the optical community as well, be- 
cause it is relatively easy to create in graphene an inho- 
mogeneous potential pattern with any distribution of p-n 
and n-n junctions, while designing a periodic or random 
stack of alternating positive-negative dielectric layers is 
nowadays a feasible task. Thus, graphene could provide 
analog laboratory models in order to test the optics of 
metamaterials in a controlled way. 

B. Dirac point 

The key feature of graphene, from which all its unique 
transport properties stem, is the existence of Dirac cones 
in the band structure of its energy spectrum. At first 
glance, one does not have to work hard to obtain in op- 
tics a double-conical, graphene-like dispersion law: it is 
inherent in any plane monochromatic wave propagating 
in a homogeneous medium, as the relation between its 
frequency ui and the wave number k is given by 

u 2 = c 2 k 2 . (1) 

However, the contact of two cones in Fig. [TJ only looks 
like a Dirac point (DP). In fact, of the two cones in Fig.[TJ 
only one (for example, the upper one) is related to a pho- 
ton, while the second solution of Eq. (UJ (lower cone in 
Fig. [IJ is redundant. This lower cone does not carry any 
additional information and is not related to any differ- 
ent physical entity like, for example, a hole in graphene 
or a positron in the case of relativistic QED. In other 
words, the photon and antiphoton are identical 3 -. Hence, 
the challenge in optics is to create a structure with a 
real DP in its spectrum, so that different types of waves 
would correspond to two different cones (a sort of op- 
tical "particlc-antiparticlc" pair). Appropriate for this 
purpose are photonic crystals, in which two modes de- 
generated in a homogeneous space become split by the 
periodicityi 5 -. 

The analytical and numerical study of two-dimensional 
periodic structures (infinite rods embedded in a back- 
ground medium with a different dielectric constant) 
was carried out as early as in 1991, for square^ and 
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FIG. 1. (color online) Surface u)(k x , k y ) described by the dis- 
persion equation (1) for electromagnetic waves. The contact 
of two cones looks like a Dirac point in graphene; however, in 
contrast to graphene, for light both cones correspond to the 
same field. 

triangular— lattices. Linear singularities (that nowa- 
days are called Dirac points) are clearly seen in the band 
structures of both systems, although they escaped the 
attention of references [l6| and [l?] mostly interested in 
absolute band gaps for different polarizations. After- 
wards, during more than two decades, studies of two- 
dimensional photonic crystals were primarily aimed on 
maximizing the photonic band gapiSr— (see also the re- 
view [2l| and references therein), until the discovery of the 
unusual transport properties in graphene, and the poten- 
tial to reproduce them in optics switched efforts towards 
the search and further exploration of photonic structures 
with Dirac-cone-likc singularities in their transmission 
spectr a 5 ' 7 ' 15 ' 22 . A number of new optical phenomena 
arising due to the existence of Dirac points were pre- 
dicted and discovered: diffusion- like 1/L-dependence of 
the pulse intensity on the distance L of propagation in- 
side the photonic crysta l 15 ' 23 (which is unusual for non- 
random media); oscillatory motion of a Gaussian beam 
(optical analog of the Zitterbewegung effect )2i^2S; extinc- 
tion of coherent backscatteringS 7 -; conical diffraction 2 ^; 
as well as the existence of graphene-like and novel edge 
states 2 ^. 

From the above-mentioned publications one can con- 
clude that the existence of Dirac points in two- 
dimensional periodic structures is a rather universal 
phenomenon, in the sense that they appear irrespec- 
tive of sample details, such as the shape and dielectric 
parameters of the "atoms" and its structural symme- 
try. For example, in the band structures presented in 
Refs. 0, EH, E3, an d HH, DPs show up at square, trian- 
gular, and honeycomb lattices. The general criteria for 
the existence of DP in periodic dielectric samples were 
discussed in Ref. [H. 

The situation in layered periodic media is quite differ- 
ent and less studied. As we show bellow, DPs cannot 
exist in periodically-layered dielectric structures built of 
mono- type (i.e., with either all positive or negative re- 
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fractive indices) dielectrics, no matter their period, size, 
and dielectric contrast between the layers. In Refs. [32| 
andHH a one-dimensional periodic array of metallic unit 
cells of a special shape was considered, which exhibited 
Dirac points created by the accidental degeneracy of two 
modes. To create photonic Dirac cones in one dimension, 
both normal and metamaterials should be used. Inter- 
estingly enough, a single DP can exist in a homogeneous 
dispersive meta-medium at a frequency at which both the 
dielectric permittivity and the magnetic permeability ap- 
proach zero simultaneously^. Eigenwaves in structured, 
infinite periodic systems built of alternating normal and 
left-handed dielectric layers, and the transmission and 
reflection from finite samples, were analyzed in Refs. l34l - 
l36l . It was shown that, when the spatial average of the 
refractive index over the period was zero, the band struc- 
ture consisted of gaps at all frequencies except for a set 
of isolated points (discrete modes) , for which the optical 
thicknesses of the adjacent layers were equal to the same 
integer number of half-wavelengths. 



C. Brief summary 

Here we study the transport properties of layered 
periodic dielectric systems "electronically similar" to 
graphene, in the sense that they possess Dirac cones 
in their photonic band gap structures. In Sec. II, we 
demonstrate, using a simple example, and discuss the 
similarity and differences between the Maxwell and Dirac 
equations as well as between the corresponding bound- 
ary conditions. Section III presents the transmission 
through potential barriers created in graphene by apply- 
ing step-like electrostatic potentials, in comparison with 
light propagation in slabs of normal dielectrics and meta- 
materials. In Sec. IV, the photonic band-gap structures 
of periodically-layered dielectrics are compared with the 
structure of the electron energy zones of graphene sub- 
jected to a periodic potential. It is shown that, unlike 
two-dimensional photonic crystal, Dirac cones in a lay- 
ered periodic medium can exist only when the medium 
consists of alternating slabs of left- and right-handed di- 
electrics (mixed samples). In Sec. V, we study the prop- 
agation of Gaussian beams of light through periodically- 
layered mixed samples built of alternating slabs with pos- 
itive and negative refractive indices, and of beams of 
charge carriers through finite graphene superlattices. We 
demonstrate the anomalous, diffusion-like dependence of 
the intensity on the propagation distance, and an analog 
of the Zitterbewcgung effect in a wide range of parame- 
ters. New unusual transport properties of such samples 
are predicted. In particular, it is shown that two tangent 
Dirac cones manifest themselves differently: given two 
beams with frequencies belonging to different cones, one 
is focused and another is defocused. The magnitude of 
the shift of the focus is independent on the distance of 
the sample from the focal plane of the incident beam and 
is proportional to the width of the sample. At oblique 



incidence, a Gaussian beam is not displaced along the 
sample even at non-zero values of the mean value of the 
dielectric constant. 



II. EQUATIONS AND BOUNDARY 
CONDITIONS 



The dynamics of the charge carriers in an external po- 
tential u in graphene is described by a spinor 

V> = (lpA,1pB) T 

whose components are related to two sublattices in the 
unit cell of the crystal^. In the low-energy limit, near 
the Dirac point, the components of this spinor obey the 
Dirac equations. When the energy w of the charge carrier 
is fixed, then the time dependence of the spinor is given 
by exp(—iwt/H), and these equations can be written as 
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w — u(x, y) \ dx 
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u(x,y) \ dx 



dy 

■ dip a 

' dy 



(2) 



where u(x,y) is the electrostatic potential, and vp is the 
Fermi velocity. 

In order to link this to Maxwell equations, we now con- 
sider, as an example, a TE electromagnetic wave (where 
the magnetic field has only one non-zero z-component) 
propagating in a homogeneous medium, and introduce 
two complex- valued functions 



£ — E y — iE x 



U = ZH Z 



(3) 



where Z — fi/e is the medium impedance, fi and e are 
the medium permeability and permittivity, accordingly. 
In Ref. [3, instead of Eqs. ©, the Reimann-Silberstain 
vector wave functions were used to derive a quantum- 
mechanical, matrix form of the classical wave equations 
in the general case of arbitrary electromagnetic fields 
propagating in media with space-dependent permittiv- 
ity and permeability. The analogy with the relativistic 
Dirac equations was noted. 

It is easy to show that for monochromatic fields £ 
and H [the time dependence is given by exp(— iurt)], the 
Maxwell equations yield 



H 



£ = - 



d£ 

kon \dx 

an 



kon 



_ M 
dy 

i— 

dx dy 



(4) 



Here, n is the medium refractive index, and fco = uj/c. It 
is evident that after the replacement 

£ <-> tpA, 

U <-> 1p B , 

nco O (w — u)/h, 

c <R- vp , (5) 
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equations (jlj coincide with Eqs. ([2]). Namely, the 
2D Maxwell equations for the complex effective fields 
Eqs. ([3]) in a homogeneous medium, and Dirac equations 
for the wave functions of the charge carriers in graphene 
become identical. The role of the refractive index of the 
corresponding effective medium is played by the quan- 
tity n e s = (w — u)/huj. Therefore if, for example, the 
potential is a piecewise-constant function of one coordi- 
nate, the corresponding graphene supcrlattice models a 
layered dielectric structure 1 ^. In particular, a layer, in 
which the potential u exceeds the energy w of the parti- 
cle, w — u < 0, is similar to a slab with negative refrac- 
tive index n. This means that a junction of two regions 
having opposite signs of w — u is similar to an interface 
between left- and right-handed dielectric media with the 
refractive indices ri\ and 712, if 



Hi 

"2 



Ui 



u 2 



(6) 



Because of this similarity, a p-n junction can focus Dirac 
electrons in graphene 3 ^ in the same way as the focus- 
ing of electromagnetic waves by the boundary between 
a normal dielectric and a metamateriali^. However it 
is important to realize that, as it follows from Eq. ((5]), a 
change of w necessarily implies the corresponding change 
of the ratio 774/712. Thus, to model the same (i.e., with 
fixed values of u% and 112) graphene bi- layer structure, 
but at different energies, one has to chose different pairs 
of dielectrics. 

It follows from Eq. © that the energy spectrum of the 
charge carriers in graphene in a homogeneous potential 
u = const is linear near k = 0, i.e., consists of two cones 
tangent at the Dirac point (Fig 1): 



(7) 
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After substituting (JS), Eq. j7]) looks exactly like the dis- 
persion law ([1]) for photons. Whilst Eqs. (2) and (4) are 
akin, the similarity between the two problems is not com- 
plete. First, unlike Dirac wavefunctions, the genuine elec- 
tromagnetic fields E and H are real, i.e., equal to their 
complex conjugates. Due to this, positive and negative 
frequencies (upper and lower cones in Fig. [TJ correspond 
to the same fields, in contrast to a Dirac spinor, which 
describes electrons, when (w — v) > 0, and holes, when 
(w — it) < 0. One further distinction is that the electric 
field E satisfies the continuity condition 



dbiE 



dE x , dE v 



da 



dy 



0, 



(8) 



which is not required for the Dirac wave functions. 

Essentially different are also the boundary conditions 
for the Dirac wave functions at the interface between two 
half-spaces with potentials U\ and u 2 , 



i'lA = ip: 



'2 A, 







IB 



0: 



2B, 



(9) 



and for the effective electromagnetic fields £ and H at 
the boundary between two dielectrics: 



Z2H1 = Z1H2, Re£i = Re£ 2 - 



(10) 



Equations (TTU)) follow from the boundary conditions for 
the real fields, E\ v = E2 V , H\ z = H2 Z (the boundary is 
assumed to be parallel to the y-axis ). 

Although the charge transport in graphene and the 
light propagation in dielectrics are governed by similar 
equations and Q inside the medium, the bound- 
ary conditions (|9|) and (|10[) arc distinct from each other. 
This means that, generally speaking, the coupling be- 
tween adjacent samples in these two cases is differ- 
ent. This difference is more conspicuous in the partic- 
ular case of a monochromatic plane wave, % cx £ cx 
exp [i {k x x + k y y — Lot)], when the boundary conditions 
for the effective fields take the form: 



£ 1 [1 + iky/k 2x ] 



Z2W1 — Z1H2, 
£2 [1 + iky/kix] 



(11) 



where k ; 



^jx — y k^n 2 - — ky. It is easy to see that Eqs. © 

and pip coincide only when k y = (normal incidence) 
and Z\ = Z2- That is, in this particular instance, the 
transmission of Dirac electrons through a junction is 
similar to the transmission of light through an interface 
between two media with different refractive indices but 
equal impedances. In other words, at normal incidence, 
any junction in graphene, either n-n, p-p, or p-n, is anal- 
ogous to a contact between two perfectly-matched di- 
electrics or microwave elements. Such an interface is ab- 
solutely transparent to the normally-incident radiation 
and therefore to the Dirac electrons in graphene as well. 
This provides a more intuitive insight into the physics of 
the Klein paradox (perfect transmission through a high 
potential barrie r 37 ' 40 ) in graphene systems. 



III. TRANSMISSION THROUGH POTENTIAL 
BARRIERS AND DIELECTRIC SLABS: 
SIMILARITIES AND DIFFERENCES 

To compare the transport properties of Dirac electrons 
and light in the general case of oblique incidence, k y =/= 0, 
we first consider the transmission of particles through a 
step-like potential barrier, i.e. through the line x = 
separating two domains (1 and 2) of a graphene sheet 
with different values u\ and «2 of the potential. In what 
follows, we will consider only one spinor component, say 
-0,4 =ip, because ipB could be found using Eq. @. The 
solutions of Eqs. ^ in both domains can be presented 
as a linear combination of plane waves with equal (at the 
chosen geometry of the system) values of k y : 



i\k y y-{w-Uj)tm 



(+) f ,ik j 
3 C 



.7 = 1,2. (12) 
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Prom now on we consider the range of parameters where 
there is no total internal reflections and therefore Imk x = 
0. From the continuity of the wave functions at the 
boundary, Eq. Q , it follows (the incident wave now prop- 
agates from medium 1 to medium 2): 



4 +) 
4~ ] 



M 



l->2 



n 
4 



(+) 



where the transfer matrix M\^ 2 for the interface x 
is equal to 



Mi_ 



1 



2 cos( 



9 tx 



9& 



and the matrix elements are 

- i02 ± Sl s 2 e ±i01 . 



o (±) 



(14) 



(15) 



Here Sj = sgn(iu — Uj), 9\ and 8 2 are the angles 
of incidence and refraction, respectively, while 0j = 
arctan(/Cy J kjx ) • 

We determine the transmission T and reflection R coef- 
ficients as the ratios of the normal-to-the-boundary com- 
ponents of the densities of the transmitted, J 2 , and rc- 
fleeted, j£ , currents divided by the normal component 
of the incident current density : 



1 ~ J 2x l J lx i JX ~ J lx l J lx ' 



where 



^=±2^008 0,- 







(±i 



(16) 



(17) 



From Eqs. (fTo]) , p7|) the following formulas can be ob- 
tained (see, i.e. ) 41 ' 42 ): 

2 cos Q\ cos 02 



T 



R = 



1 + cos (s 1 9 1 + s 2 9 2 ) 
1 - cos (si#i - s 2 9 2 ) 
1 + cos (si0i + s 2 6 2 ) ' 



(18) 



where the angles of incidence 9\ and refraction 9 2 are 
connected by the relation 



sm t/9 = sm 9\ . 

W — 1t 2 



(19) 



The same signs of s\ and s 2 correspond to n-n (or p-p) 
junctions, while for a n-p (or p-n) junction s\ = —s 2 . It 
is easy to see that at normal incidence {6\ = 9 2 = 0) the 
potential barrier of any height is absolutely transparent 
at any energy (the Klein tunneling effect). 

In the case of light propagating through the interface 
between two dielectrics, whose parameters are £i, /ii and 
e 2 , (j,2i the transfer matrix 



Mi. 



l 



2 cos 9 2 



q{+) 



{-) 



»l->2 
g{+) * 
»l->2 



(20) 



connects the amplitudes of leftward and rightward prop- 
agating waves at both sides, and has the same form as 
Eq. (fT4")l , with the matrix elements g 2 ± ^ 1 replaced by 



Q\4i = cos 



9 2 ± s 1 s 2 cos9 1 ^- 



(13) 

for TE waves, and by 



yi_/ 2 = cos6>2 — ± sis 2 cos 9 1 
"I 



(21) 



(22) 



for TM radiation. In Eqs. (|2T|) 
±yjEjfjLj, Zj = y/fij/ej. For left-handed dielectrics (e < 
0, /i < 0), the refractive index n is negative. 

When Z\ = Z 2 , the expressions (|2 1 1) and (|22|) are iden- 
tical and are related to the matrix elements (fT5")) of the 
corresponding transfer matrix in graphene, Eq. (|13[) . as 



sgnn 



q{±) 
»l->2 



Reg[% 



(23) 



This relation is a consequence of the abovementioncd dif- 
ference between the solutions of the Dirac and Maxwell 
equations: the former are complex-valued functions, 
while the electromagnetic fields are real. 

The light transmission and reflection coefficients are 
determined as the ratios of the normal components of the 
transmitted (for T) and reflected (for R) energy fluxes 
divided by the normal component of the incident energy 
flux. When Z\ = Z 2 (the situation most favorable for 
the analogy between Dirac electrons and light) T and R 
take the forms: 

4 cos 8\ cos 9 2 
= ~ : ~T2 



R = 



(cos 01 + COS0 2 ) 
(cos 01 — COS0 2 ) 
(cos 01 + COS0 2 ) 



(24) 



One can see that even in the particular case, Z\ = Z 2 , 
the equations (fTg)) and (f2"4"]l are, generally speaking, dif- 
ferent, and coincide only when 9\ = 9 2 = 0, i.e., when 
the boundary conditions Eq. pip arc equivalent. This 
means that in spite of the identity of the equations @ 
and ((H) , the analogy between the transport of Dirac elec- 
trons in graphene and electromagnetic radiation in di- 
electrics should not be extended too far. Due to the dif- 
ferences in the boundary conditions, the analogy holds 
only for normal incidence on the interface between two 
perfectly-matched media. 

The transmission of Dirac electrons trough a potential 
barrier of finite width d, and of an electromagnetic wave 
through a dielectric slab of the same width are described 
by similar matrices of the form 



B = A 2 



iSdAi 



(25) 



where the indices 1 and 2 now correspond, respectively, 
to the outside and inside of the barrier (slab). The 
matrix A is equal to M, as in Eq. (fl4"|) for graphene, 
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and A = A4, as in Eq. ((20)) for dielectrics. Obviously, 
yl.2-s.i^4i->.2 = I, where I is a unit matrix. The diagonal 
matrix Sd = diag (e lv , e~ lip ), with ip = k 2x d describes 
the propagation inside the barrier (dielectric slab). From 
Eq. (|2"5)l it follows that under the condition 



where 



k 2x d ■ 



m = 1, 2, 3, 



(26) 



B = (— l) m I. This means that either a potential bar- 
rier or a layer of dielectric is transparent if its width is 
equal to an integer number of half-wavelengths. This 
is a general property of both, Dirac and Maxwell equa- 
tions, which is independent of the ratio between the 
impedances, and its physical nature has nothing to do 
with the Klein tunneling. 



IV. PERIODICALLY-LAYERED GRAPHENE 
SUPERLATICES AND PHOTONIC 
STRUCTURES 

In this Section, we compare the structure of the elec- 
tron energy zones of graphene subject to a periodic elec- 
trostatic potential, with the photonic band gap struc- 
tures of periodically-layered dielectrics. To do this we 
assume that the potential in graphene takes the two val- 
ues u\ and u 2 in alternating areas of the thicknesses 
d\ and d 2 , and the corresponding dielectric sample is 
built of alternating layers of the same thicknesses, d\ and 
d 2 , with refractive indices n\ and n 2 , and impedances 
Z\ and Z 2 , respectively. In both cases, the propaga- 
tion in the layers 1 and 2 is described by the matrices 



diag (e^i , e" 



1,2, where ipj = djkj 



Assuming that the layers are parallel to the y-axis, the 
transformation matrix P on the period D = d\ + d 2 is 
defined by tjj(x + D) = P^{x), and is equal to 



P = S 1 M 2 ^ 1 S 2 M 1 ^ 2 . 



for graphene, and 



2^1 S 2 M l->2- 



(27) 



(28) 



for dielectrics. 

The eigenvalues A = exp(ik\\D) (k\\ is the Bloch 
wavenumber along the x-axis) of the matrix P (V) de- 
pend on the energy (frequency) and on the tangent com- 
ponent k y of the wavevector. Both periodic structures 
are transparent if 



1. 



(29) 



In other words, Eq. (|29[) determines the transparency 
zones: the ranges of the energies (frequencies) and wave 
numbers k y for which the longitudinal wavenumber fc|| is 
real^. It can be shown that 



A = F ± \/.F 2 - 1, 



(30) 



tan B\ tan 9 2 — 



F(w, k y ) = cos ip i cos tp 2 
sis 2 



shupi sin tp 2 (31) 



cos f i cos V 2 / 
for Dirac electrons in graphene and 

F{u>, ky) = cos ipi cos ip 2 

sxs 2 /Zicosfli Z 2 cos6> 2 \ . , . 

tz 7, — r -n sin w\ sm w 2 . 32 

2 \Z 2 cos6 2 Zicosej * y V ' 

for electromagnetic waves in layered dielectrics. In both 
cases, the eigenmodes obey the dispersion equation: 



cos(/c||£>) = F. 



(33) 



It follows from Eqs. (|29|) and ([30]) that a periodic struc- 
ture is transparent for the points in the plane {k y ,w) [or 
(k y ,uS)\ for which the inequality |F| < 1 holds. One the 
other hand, in each of these two planes, the conditions 



(px=pn, p= 1,2,3,. 
ip 2 = qir, q = 1,2,3, 



(34) 



(jj and q are positive integer numbers) determine two sets 
of curves, where P and V arc equal to (— l) p+q I. This 
means that those types of curves belong to transparency 
zones. It is apparent that since at the crossings of these 
lines, the eigenvalues of the matrices P and V are equal 
to ±1, such singular crossing points lie at the edges of 
the zones and represent their points of contact, known as 
"diabolic points" (see, i.e., Ref. HH). 

The coordinates, k v t and fcot = ^t/c, of these points in 
the plane {k y , kg) can be found from Eqs. (|34|) : 



k 2 



(pn/di) 2 - (qir/d 2 f / (n\ - nf) 



1-2 

K yt 



n\ {p-n/dif - n\ (qir / ' d 2 f / (n\ - n%) . (35) 



In the case of graphene, the coordinates of the diabolic 
points in the plane (w, k y ) are calculated in the same way. 
In the vicinity of the points Eq. (|34p . the phases ipi and 
ip 2 can be written as 



<pi = pn + S(fi 



where \5tfj\ <C 1. Substituting Eqs 
yields: 



f 2 = qn + 5ip 2 , (36) 
into Eq. 



1 



1 — — (5ip\ + Stpl + 2a8ipi8ip 2 ) 



1 - 2^ S(fl,Slf2 ' > 



(37) 



with 



SiS 2 ( ZiCOS^i Z 2 cos 9-2 



Z 2 cos 9 2 Zi cos 6*i 
_ 1 / ki x e 2 k 2x e 1 
2 \k 2x £i k\ x e 2 
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where the components k\ x and k% x are taken at the point 
(k yt ,kot). In the plane (Sk y ,Sk ), where 5k and 6k y are 
small deviations of kg and k y from their values given 
by Eq. Q35p . the points for which the quadratic form 
f(S(pi,S(p2) P7[) is positive, constitute a transparency 
zone, while the points with / < correspond to a gap in 
the photonic spectrum. 

The relation between Sipj, 5k y , and Sko follows from 
the formula: 

<Pj=djyJn}(ko)l<%-k*, (38) 

where nj{k$) is the refraction index of the j-th layer. An 
analogue formula for graphene has the form: 

fj = d 3\J[{ w ~ u j)l hv F] 2 ~kl, (39) 

While in conventional dielectrics, the dispersion can 
be ignored, if it is weak enough, in left-handed ma- 
terials it always must be taken into account. Indeed, 

the surface oj(k), where k = ^Jk x + k y , for normal di- 
electrics, is a cone similar to the one presented in Fig. [T] 
i.e., uj(k) = n~ 1 ck, with n = const > 0. For left- 
handed metamaterials n < 0, and the group velocity 
v g = (diu/dk) is negative, v g < 0, i.e., is antiparallel 
to the phase velocity tu/k. The surface u(k), in a small 
vicinity of a certain frequency ujq, has the form depicted 
in Fig. [5] Although locally it looks like a part of the 
lower cone in Fig. [1] uj(k) ^ n~ 1 ck, with n = const on 
this surface, because the tangent point of the cones is 
shifted from the origin. 




FIG. 2. (color online) Surface u)(k x , k y ) in a small vicinity of 
a certain frequency for left-handed media. 

Although the absolute values of the transmission and 
reflection coefficients in Eqs. (jT5J) and (|24l) are different, 



the dispersion characteristics of two adjoining right- and 
left-handed dielectric layers, and the energetic spectrum 
diagram of n-p junction in graphene are identical, as it is 
shown schematically in Fig. [3] As a consequence of this 
identity, the periodic dielectric structure and graphene 
superlatticc formed by a periodic external potential pos- 
sess the same unique transport properties. 
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FIG. 3. (color online) The dispersion characteristics cj(fc), 
shown as blue lines, of two adjoined right-handed (left panel) 
and left-handed (right panel) dielectric layers is identical to 
the energetic spectrum w(k) of n-p junction in graphene. The 
horizontal dashed red line represents the wave frequency (left 
panel), or the quasi-particle energy (right panel). 



A. Photonic structure 
From Eq. ([38]) we obtain: 

S(pj = f^lSko - k yt Sk y ) , (40) 

kjx V v jg / 

where Vj g is the wave group velocity in the j-th layer 
at the frequency u> t = c&ot- The group velocity in (|40[) 
is positive for normal dielectric layers, and negative for 
layers of left-handed metamaterials. 

Equation (|40|) shows that the phase variations Sipj van- 
ish on the lines 

Sk = l ^Sk y . (41) 
ck Jt 

These lines lie in the transparency zones and intersect at 
the point of their tangency. The line associated with the 
right-handed dielectric lies in the I-st and the Ill-rd quad- 
rants, whereas the line associated with the left-handed 
metamaterial lies in the Il-nd and the IV-th quadrants, 
as it is shown in Fig. |H When the lattice is composed of 
conventional dielectrics only, both lines 5(fi(5k y , 5ko) = 
and Sip2(Sk y ,Sk ) = lie in the I-st and the Ill-rd quad- 
rants, and the transparency zones touch one another in 
a manner shown in Fig. 3^.. There are band gaps above 
and below the frequency ui t = ckat, and these gaps are 
degenerate in point-like non-transparent zones (gaps) at 



8 



the frequency u = u>t- Note that the zone structure of the 
lattice composed of both left-handed dielectrics is similar 
to the one shown in Fig. 0^,, but symmetrically reflected 
with respect to the 5ko axis. 






SK (b) 
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FIG. 4. (color online) Fransparency zones (gray) in the 
plane oj(k x ,k y ) for: (a) a lattice composed of conventional 
dielectrics only; (b) and (c) lattices containing both left- and 
right-handed dielectrics. In (b), the point of the zones tan- 
gency presents a point-like transparent zone. In (c), the con- 
tact of two zones is a point-like non-transparent gap. Fhc 
phase deviations 6ipi z 2 vanish on the red dashed lines. 

The band structure can be significantly different when 
the lattice contains both left- and right-handed di- 
electrics. In this case, there are transparency zones above 
and below the frequency u> t = cfcot, that touch one an- 
other at the frequency uj = Wt, forming a point-like trans- 



parent zone as it is shown in Fig.|4}D. This band structure 
is similar to the energy spectrum of the charge carriers in 
graphene, and manifests genuine Dirac points. It must 
be emphasized that for such points to exist, the layers 
with both positive and negative refractive indices must 
be present, the average (over the period) value of n should 
not necessarily be equal to zero. Materials with effective 
e and/or effective fi near zero have become recently the 
subject of intensive investigation due to their unusual 
transport properties, including the existence of conical 
singularities in the band gap structures^ - — , however this 
topic lies outside the domain of our paper. 

It is worth noting that the zone structures shown in 
Fig. @] have been calculated for TE (p-polarized) waves. 
In the case of s-polarizcd fields, the dielectric permit- 
tivities £j in Eq. (|37[) arc replaced by the magnetic per- 
meabilities pLi, therefore the photonic band structure is 
different, while the coordinates of the diabolic points are 
the same. This means that for a given frequency and 
angle of incidence, the same sample could be transpar- 
ent for s-polarization and opaque p-polarization, and vice 
versa. This property can by utilized for the separation of 
polarizations. 

The transparency zones depicted in Figs. [4] are the 
projections onto the plane (k y ,ko) of the surface u = 
w(fc x , k y ), described by the dispersion equation and 
shown in Fig. [5] The zone structures, presented in 
Figs. |4^,,c and [4}d, correspond to the surfaces shown in 
Figs.[Sk and[SjD, respectively. 




FIG. 5. (color online) Surfaces 5ko(dk», Sky) for: (a) the zone 
structure depicted in Fig. [4b; (b) the zone structures shown 
in Figs. 2K,c. 

It is important to note that the presence of both left- 
and right-handed dielectrics in the periodic structure is a 
necessary, but not sufficient condition for the graphene- 
like band structure to exist. Generally speaking, the lines 
(l4"TT) and the band edges do not coincide, and can be 
widely separated. When the distance between the lines 
(1411) and the zone edges is large enough, as, for example, 
in Fig. St, the zones shape is similar to that of a lattice 
formed by right-handed dielectrics only (see Fig. 0J;) . 

The fundamental qualitative difference between the 
zone shapes of mono- and mix-lattices owes its origin to 
the strong dispersion, which is an inherent characteristic 
of left- haded media. Ignoring this fact ( "for simplicity" , 
as this is sometimes done) leads to wrong results for the 
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zones structure in the vicinity of their tangency. 



B. Graphene superlattice 



It follows from Eq. ([39]) that in a graphene superlattice 
created by an electrostatic potential with periodically- 
alternating values U\ and 1*2, the phase variations 8tp 
vanish on the lines 



h 2 v 2 k . 
w t - u. 



(42) 



where Wt and k yt are the coordinates of the zone tan- 
gency point in the plane (k y , w). The slopes of these 
lines have opposite signs when the energy w t lies be- 
tween the values u\ and 1*2 of the potential, i.e., when 
sgn[(w — ui)(w — u 2 )] = — 1. In this case, the zones touch 
each other either as it is shown in Fig. |4p or Fig. |4j: (de- 
pending on the relation between the layer thicknesses d\ 
and d,2 1 and the tangency point indices p and q). In the 
opposite case, when sgn[(w — ui)(w — 1/2)] = +1, the 
zones structure is similar to the one shown in Fig. [4^,, 
or symmetrically reflected with respect to the ordinate 
axis. Thus, the point-like transparent zones (new Dirac 
points) can appear in the graphene superlattice only in 
the energy range between u\ and u 2 - 



TRANSMISSION NEAR THE DIRAC 
POINTS 



As it was mentioned in the Introduction, the similar- 
ity between the energy spectra of electromagnetic waves 
in homogeneous media, and Dirac quasiparticles is only 
formal and physically meaningless, because the two cones 
in the photon spectra are identical. To demonstrate that 
the singular points (presented above) in the band gap 
structure of mixed periodic dielectric samples possess the 
properties, which make them true Dirac points, we con- 
sidered the transmission of a monochromatic wave of a 
frequency to through a finite stack of alternating left- 
and right-handed dielectric slabs. The dependences of 
the amplitudes and phases of the complex transmission 
coefficients t{k v ) of the y-component of the wavevector 
[i.e., of the angle of incidence 9 = arcsin(cfc y /o;)] are 
shown in Figs. |6^i and |6p, for two different frequencies 
belonging, respectively, to the upper and the lower cones 
in Fig. [5] While the amplitudes |t(A: y )| are similar, the 
phases (3 = &rg t(k y ) manifest quite distinct behaviors: 
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FIG. 6. (color online) The magnitude (blue) and the phase 
(red) of the complex transmission coefficient t(k y ) of a 
monochromatic wave propagating through: (a, b) finite pe- 
riodic stacks of alternating left- and right-handed dielectric 
layers; (c) a stack of normal dielectric layers. The frequencies 
in (a) and (b) belong to the upper and lower cones in Fig. [5] 
respectively. The k y component of the wavevector determines 
the angle of incidence 9 = arcsin (ck y /ui). The parameters of 
the numerical simulations used here are: number of periods 
AT = 9, p = 1, q = 8, di = 0.35D, ei = m = 1, |e 2 | = 0.8, 
|n 2 | = 2.5, \v 2a \ = 0.25c, and D\Sk \ = 0.05. 



j3{hy) — lb6 {hy kyjy^) 



(43) 



linear) function of the angle of incidence at all frequen- 



where opposite signs correspond to two different cones, cics, and has no singularities at k y 



and k ym is the position of the phase fi{k y ) extremum. 
For comparison, the dependences of and /3(k y ) for 

a periodic stack of normal dielectric layers are shown in 
Fig. IHt- In this case, the phase is a monotonic (close to 



Exactly the same parabolic dependence of the phase 
/3 (k y ) of the transmission coefficient is intrinsic to 
graphene superlatticcs, when the point of tangency cor- 
responds to the point-like transparent zone. The de- 
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pendences of the amplitudes and phases of the complex 
transmission coefficients t(k y ) of the y-component of the 
wavevector are shown in Fig. [7] for two different energies: 
one above (Fig. [7^,) and another below (Fig. [TJd) the en- 
ergy wt of the zones tangency. 
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FIG. 7. (color online) Magnitude (blue) and phase (red) of 
the complex transmission coefficient t(k y ) of a monoenergetic 
quasiparticle propagating through: (a, b) finite periodic set of 
alternating p-n and n-p junctions in graphene. The quasipar- 
ticle energies are above (a) and below (b) of the energy of the 
zones tangency Wt- The indexes of the zones tangency point 
are p = 1, q = 5, and the layer thicknesses are di = 0.3D and 
d,2 = 0.7D. In the same lattice, the zones tangency point, 
whose indexes are p = 3 and q = 1, corresponds to a point- 
like gap, and the dependence t(k y ) looks like the one for a 
periodic stack of normal dielectrics (Fig. [BJi) 

As it was mentioned above, the presence of both left- 



and right-handed dielectrics, or alternating p-n and n- 
p junctions in graphene, is not sufficient for a point of 
zone tangency to be a point-like transparency zone (Dirac 
point). In the same lattice, points of the zone-tangency 
with different indexes p and q can be either point-like 
gaps, as shown in Fig.[7J;, or point-like transparent zones, 
i.e., Dirac points (Figs. [7^,b). 

A number of interesting effects emerge in the vicinity of 
the Dirac point of the mixed periodic sample or graphene 
superlattice. Some of them are caused by the parabolic 
dependence (|4"5)) of the transmittance coefficient phase of 
the angle of incidence (of k y ). Let us consider the prop- 
agation of a monochromatic beam of light, bounded in 
the transverse dimension (Gaussian beam, for instance) , 
through the mixed stack of a thickness L. Generally, a 
beam transmitted through a slab of a normal dielectric 
is shifted along the surface, as it is schematically shown 
in Fig. [5^. The shift A is determined by the phase /3 of 




FIG. 8. (color online) (a) A beam of light transmitted through 
a slab of a normal dielectric is shifted along the surface. The 
shift of the beam passing through a mixed periodic stack is 
much smaller, can be either positive or negative, or even equal 
to zero, as shown in figure (b). 



the transmission coefficient t: 

A _ dP{ky) 



(44) 



™ y — "i/O 



where k y o is the tangential component of the wave vector 
of the central ray in the incident beam. It is assumed 
in Eq. (|4"4")) that the angular width Ak y of the incident 
beam is rather small. Since mixed periodic stacks are 



characterized by the parabolic dependence /3(k y ), (|43[) . 
d/3(k y )/dk y is small or even equal to zero, in which case 
the longitudinal shift A is absent completely. 

The absence of the longitudinal shift, when 
dj3{k y ) / dk y \ k _ k a = 0, can also be observed in 

the so-called near-zero-index metamaterials^ and in 
ID periodic lattices with birefringent materials^. It is 
important to note that in our case, this phenomenon has 
a different physical origin. 

The parabolic dependence f3(k y ) affects also the cur- 
vature of the transmitted beam phase front, i.e., it shifts 
the focal plane of the incident Gaussian beam. The value 
A f of this shift is proportional to the thickness L of the 
mixed stack and is independent of the stack position on 
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the beam trajectory. The sign of the shift Af depends 
on the sign of the detuning of the beam frequency w from 
the Dirac-point frequency u) t . The focal plane is shifted 
forward (defocusing) when ui < uit, and backward (focus- 
ing) when uj > oj t . 

This surprising focusing properties of mixed periodic 
samples are demonstrated in Fig. [§] Of the two beams 
with the frequencies resting on two different dispersive 
cones, the one corresponding to the upper cone is focused 
by the sample, blue curve, while the other (lower cone) 
is defocused, red line. For comparison, the black curve 
presents the intensity distribution in the beam propa- 
gated in free space. This phenomenon is highly unusual 
by itself, and also reinforces the similarity of the tan- 
gent points of the cones to a genuine "optical" Dirac 
point: different cones are not identical, and represent ob- 
jects with distinct physical properties - a sort of optical 
particlc-antiparticle pair. 
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FIG. 9. (color online) Focusing properties of a mixed peri- 
odic sample. After the transmission through a mixed peri- 
odic sample (grey area), the focus of the beam with a fre- 
quency belonging to the upper cone (blue curve) is shifted to 
the left with respect to the beam propagating in free space 
(black curve). The same sample shifts the focus of the beam 
with frequency from the lower cone (red curve) in the oppo- 
site direction. The parameters of the numerical simulations 
used here are: number of periods N = 100, p = 2, q — 5, 
di = 0.4D, ei = m — 1, £2 = —0.5, ni = —0.5, V2 g = —0.6c, 
and D\Sk \ = 0.3. 



When the frequency oj of the incident beam belongs to 
a point-like gap of the corresponding infinite structure, 
i.e, oj = Wt, one would expect an exponential decay of 
the transmitted beam intensity Ttr as a function of the 
mixed stack thickness L: Itr(L) oc exp(— jL). However, 
the well-pronounced constant asymptotic of the function 
L-I tl -(L) (Fig. [TD1 red line), demonstrates an anomalously 
high intensity, decaying as 1/L. Such a diffusion- like de- 



pendence is one of the consequences of the linear, Dirac 
cone-like dispersion. In two-dimensional photonic crys- 
tals with triangular lattices of normal dielectric rods, this 
phenomenon was predicted in Ref. [Tol In layered struc- 
tures, it could take place only in the presence of left- 
handed elements. 
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FIG. 10. (color online) Intensity I(N) (blue line) of the trans- 
mitted beam as a function of the distance of propagation in- 
side a mixed stack. The distance is measured in the numbers 
of periods N. The frequency u of the incident beam belongs 
to a point-like transparency zone (Dirac point) of the corre- 
sponding infinite structure, uj = u) t . The large- iV asymptotic 
of the red line, NI(N), is constant, which means that the in- 
tensity is inversely proportional to the distance (diffusion-like 
dependence). The parameters of the numerical simulations 
here are: p = 1, q = 8, dx = 0.35D, ex = nx = 1, £2 = —0.8, 
?i2 = —2.5, and V2 g = —0.25c. 

Because of the similarity between the energy spectrum 
of relativistic electrons and the frequency band structure 
of a mixed periodic dielectric structure, it is natural to 
assume that the light beam propagation into mixed sam- 
ples can be accompanied by a Zitterbewegung-likc phe- 
nomenon. Indeed, the spatial distribution of the energy 
flux inside the mixed sample manifests a trembling mo- 
tion: the "center of gravity" of the flux oscillates in the 
transverse direction (along the y-axis). In Fig. [TTJ the 
spatial distribution of the energy flux inside the mixed 
sample is shown. The oscillatory motion of the center of 
gravity I c (y) is clearly seen in Fig. [12] Note that Fig. QT] 
is similar to Fig. 3 in Ref. [52l. where the probability func- 
tion of a moving electron (solution of the Dirac equation) 
is shown. 

Figure [11] also demonstrates the two above-mentioned 
effects: the absence of the longitudinal shift (the energy 
propagates normally to the sample boundary, unlike the 
angle of incidence of the beam is far from the normal) 
and the focusing of the beam. 

All these effects — the absence of the longitudinal 
shift, the focusing of the beam, and the Zitterbewegung 
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FIG. 11. (color online) Spatial distribution of the energy 
flux of the beam propagating (from left to right) through the 
mixed sample. The boundaries of the sample are marked 
by the vertical dashed lines. The "center of gravity" of 
the flux oscillates in the transverse (along the y-axis) di- 
rection. The parameters of the numerical simulations used 
here are: number of periods N = 100, p = 4, q — 10, 
di = 0.4D, Ei = m — 1, £2 = —0.5, ?i2 = —1.5, «2 S = —0.2c, 
DSk = -0.07, and D5k y = 0.5. 
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FIG. 13. (color online) Transmission of the Gaussian beam of 
charge carriers through a periodic set of p-n and n-p junctions 
in graphene near a point-like gap. The sample boundaries 
are marked by the vertical dashed lines. The parameters of 
the numerical simulations used here are: number of periods 
N = 100, p = 1, q = 5, dx = 0.3D, m = -ua = 20hv F /D, 
Sw = -Q.Zhvp/D, and D8k y = 0.2. 
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FIG. 12. Oscillations of the transverse coordinate of the "cen- 
ter of gravity" of the beam propagating in the mixed sample. 

phenomenon — are also seen in graphene super lattices. 
As an example, the current density distribution in the 
Gaussian beam of monocnergetic charge carriers propa- 
gating through a finite graphene superlattice is depicted 
in Fig. [131 

Let us now compare the wave propagation through a 
mixed dielectric slab or a graphene superlattice with the 
transmission of an electromagnetic wave through a plate 
made of a homogeneous (right- or left-handed) dielectric. 



In a small vicinity of the normal angle of incidence, the 
dependence (3{k y ) of the transmission coefficient phase 
(3 of the tangential component of the wave number k y 
(the dependence of the angle of incidence) has the same 
parabolic form as discussed above [see Eq. (|43pj. with 
the only difference that at homogeneous dielectric plates, 
k y t = by definition: 

P{ky) = ±bkl (45) 



Here the plus sign corresponds to right-handed di- 
electrics, while the minus sign corresponds to left-handed 
dielectrics. In the first case, the refraction is positive, i.e., 
the incident beam and the beam inside the plate lie on op- 
posite sides from the normal; whereas in the second case, 
the beams lie on the same side of the normal (negative 
refraction). Therefore, the shift A is positive for plates 
of a right-handed dielectric, and negative for left-handed 
dielectrics. In other words, in normal samples, k y and A 
have the same signs, while in mctamatcrials the signs are 
opposite. Amazingly, these two situations, each inher- 
ent to different kinds of homogeneous materials, can be 
observed in the same mixed periodically-layered sample, 
and graphene superlattice. Indeed, depending on which 
side of the Dirac point the frequency lies, the refraction is 
either positive (upper cone), or negative (lower cone). In 
this regard, it is more appropriate to refer to these two 
cones as a "medium - anti-medium" pair, rather than 
"particle-antiparticle" , as in homogeneous graphene. 
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VI. CONCLUSIONS 

We have shown that some of the exotic properties of 
charge transport in graphene can be reproduced studying 
the propagation of light through layered dielectric sam- 
ples. Similarities and distinctions between Maxwell and 
Dirac equations, and between the corresponding bound- 
ary conditions have been studied. Although the equa- 
tions for the real electric and magnetic fields are essen- 
tially different from those for the Dirac electrons, under 
some conditions they can be reduced to similar form. For 
example, Eq. (|4]) for the complex combinations given by 
Eq. ([3]) coincide with the Dirac equations Eq. $Z§. There- 
with, the role of the refractive index in graphene is played 
by the difference between properly normalized values of 
the Fermi energy and the external electrostatic potential 
[Eq. (O]. The boundary conditions for a Dirac quasipar- 
ticle incident on a plane separating two areas with dif- 
ferent potentials, and for a electromagnetic wave prop- 
agating through an interface between two layers of ho- 
mogeneous dielectrics are, generally speaking, different. 
They coincide only when the impedances of the layers 
are equal, and the direction of the propagation is normal 
to the boundary. It means that at normal incidence, any 
junction in graphene is analogous to a contact between 
two perfectly-matched dielectrics and therefore is abso- 
lutely transparent to normally-incident Dirac electrons. 
This provides a more intuitive insight into the physics of 
the Klein paradox in graphene. 

The analytical and numerical analysis of the photonic 
band gap structures of infinite periodically-layered sys- 
tems reveals an infinite number of the so-called "dia- 
bolic points" (singular points of contact of two trans- 
parency zones) in the (u>, 9) spectral diagrams (examples 
are shown in Fig. U) . A distinction needs to be drawn 
between two types of these singularities: point-like trans- 
parency zones, like in Figs.|4^, and[4}:; and point-like gaps 
in the spectrum, similar to the one presented in Fig. 
Although all three pictures in Fig. @] are topologically 
equivalent, the transport properties of the correspond- 
ing finite periodic stacks of layers differ drastically in 
the vicinities of these points. Waves with frequencies ly- 
ing on opposite sides of the singularities of the first type 
propagate through the samples in similar ways. In the 
same time, when two tangent spectral cones form a point- 
like gap, the electromagnetic radiation interacts with the 
same sample differently, depending to which cone its fre- 
quency belongs to. Studies of the propagation of beams 



of light show that only the diabolic points of this type 
posses the properties of genuine Dirac points. We demon- 
strate that in mono-type layered structures (i.e., in those 
built of either normal or left-handed dielectrics) just the 
diabolic point of the first type can exist, and conical, 
Dirac-type singularities appear only in mixed (with al- 
ternating left- and right-handed layers) samples. This is 
in contrast to two-dimensional media where Dirac points 
were discovered in various types of photonic crystals with 
normal dielectric elements. It is important to note that 
the physical nature of the Dirac points that we consider is 
different from that in systems with zero average value of 
the refractive index: in mixed layered structures they are 
due to the specific strong dispersion (phase and group ve- 
locities have different signs) inherent in the elements with 
negative refraction. 

Although the angular dependences of the transmission 
and reflection coefficients from a single interface in lay- 
ered dielectrics and graphene supcrlatices are different 
[compare formulas Eqs. (fT5)l and (|24p] . the spectral prop- 
erties of these two structures are conceptually identical 
and entail similar features in the light and charge trans- 
port. Considering, as examples, the transmission of the 
Gaussian monochromatic beams of light and monoener- 
getic Dirac electrons through the corresponding (dielec- 
tric or graphene) samples we predict the following Dirac- 
point-induced effects: (i) two tangent Dirac cones influ- 
ence the propagation of a beam in different ways: the 
beam is focused when the frequency (energy) belongs to 
the upper cone, and is defocused at frequencies (ener- 
gies) lying in the lower one; (ii) the transverse shift of 
the beam is anomalously small or even zero; (iii) the de- 
cay of the intensity at forbidden frequencies is diffusion- 
like; (iv) a spatial analog of the Zitterbewegung effect 
(i.e., trembling motion of the "center of gravity" of the 
energy flux) is observed in periodically-layered dielectric 
structures with non-zero average refractive indices and in 
graphene super-latices. 
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